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Abstract 

Measurements of jet suppression in Pb+Pb collisions by the ATLAS Collaboration are reported. The production 
of inclusive jet yields as a function of jet pj, collision centrality and jet size parameter R are measured and presented 
through the central-to-peripheral ratio, Rqp- Jets are found to be suppressed in central collisions relative to peripheral 
collisions by approximately a factor of two. The suppression is found to show almost no variation with jet p T , and the 
Rep is found to increase slighly with increasing R. Measurements of heavy quark energy loss are also presented using 
muons from semi-leptonic decays of heavy flavor hadrons. The maximal suppression is observed to be Rqp ~ 0.5 and 
shows no significant depedence on the muon pj. 
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1. Introduction 

Relativistic heavy ion collisions provide access to matter at the highest temperatures ever produced in the labo- 
ratory. In these collisions the produced medium consists of a system of deconfined quarks and gluons known as the 
quark-gluon plasma (QGP). Jets produced in these collisions must traverse the medium and may suffer energy loss or 
modification of their parton shower through interactions. This phenomenon, known as jet quenching, provides direct 
access to the QGP transport properties through the jet-medium interaction. One of the first proposed signatures of 
QGP formation was the energy loss of back-to-back jets traversing different in-medium path lengths [ 1 1. Experimental 
signatures of partonic energy loss were first observed at RHIC through the measurement of the suppression of single 
hadrons [2 3 4 5| and modification of di-hadron angular correlations EE) . Recently the phenomenon of jet quench- 
ing was observed directly using fully reconstructed jets at the LHC EI9). The measurement of highly asymmetric 
dijets in central Pb+Pb collisions can be interpreted as an experimental realization of Bjorken's picture, where the 
large asymmetry is attributed to dijet pairs where one jet is produced near the surface and suffers little modification 
while the balancing jet suffers significant energy loss in the medium. Furthermore, the lack of modification of the dijet 
angular correlation indicates that the quenching mechanism must allow for significant energy loss without distorting 
the jet's angle. 

The asymmetry measurements are sensitive to differential energy loss- the energy loss of one jet relative to another. 
Additional information about the quenching mechanism by considering the inclusive energy loss of jets, which is 



Email address: anger amiOcern. ch (A. Angerami for the ATLAS Collaboration) 



A. Angerami for the ATLAS Collaboration / Nuclear Physics A 00 (2012) 7-[§] 



2 



a measure of the average, absolute energy loss. As the jet production spectrum is steeply falling, the systematic 
downward shift in the jets' pi caused by the energy loss results in a reduction of the yield of jets at a given pj, 
a phenomenon referred to as jet suppression. The suppression can be quantified by comparing the jet pj spectrum 
in central and peripheral collisions. As central collisions involve a larger nuclear overlap in the collision, there is a 
geometric enhancement to the hard-scattering rates. To assess the effects of the quenching this enhancement must be 
removed, thus the per-event jet yields must be normalized by a factor, R co n, which corrects for the centrality -dependent 
geometric enhancement of parton luminosity per nuclear collision. The suppression is then evaluated in terms of a 
nuclear modification factor, 



Rcp(.Pt)\ 



1 AW 



?coB -^jet 
N evt 



(1) 



periph 



This quantity has been measured first at RHIC JU |3] |4] and later at the LHC ifTUl [TTI for single hadrons, however 
such measurements do not provide the full information about the quenching mechanism as they are only sensitive to 
part of the jet. At the LHC, measurements of fully reconstructed jets are possible due to both the larger acceptance 
of the detectors, especially calorimeters with large pseudorapidity converage, as well as the enhanced hard scattering 
rates for jets with energies significantly above the underlying event (UE) fluctuations. The clustering in jet algo- 
rithms such as the anti-£, |[T2l are controlled by a size parameter R, which corresponds to a nominal cone radius 
of R — y'A.rj 1 + Acp 2 about the jet axis. The angular pattern of medium-induced radiation can be investigated by 
studying the suppression of jets with different R values, and can thus provide quantitative constraints on jet quenching 
calculations fl3l . 

Heavy flavor jets provide additional insight into the quenching mechanism; the additional scale introduced by the 
large quark mass may force a different interplay between elastic and inelastic energy loss than that for light quarks [ 14 1 . 
However, measurements of heavy quark production at RHIC via semi-leptonic decays to electrons showed a combined 
charm and bottom suppression in Au+Au collisions comparable to that observed for inclusive hadron production lfT5l 
[161Q2)- There is disagreement in the theoretical literature regarding the interpretation of the RHIC heavy quark 
suppression measurements lfl8l [T9l l20l l2T1l regarding the role of non-perturbative effects ll22l l23l l24l . Heavy quark 
jets can be studied by considering the semi-leptonic decay process of heavy flavored hadrons within the jet. At LHC 
energies this process is the primary source of single muon production at intermediate pj values (4-14 GeV) (25). Thus 
the suppression of single muons provides an indirect measurement of the jet suppression of heavy quark jets. 

In these Proceedings the suppression of inclusive jets reconstructed with the anti-£, algorithm using R values, 
R = 0.2, 0.3, 0.4 and 0.5 and the jet Rqp is presented as a function of pj, R and centrality, and the discussion 
closely follows the results of a recent ATLAS paper ll26l . All results are fully unfolded for experimental effects 
using an unfolding procedure based on the singular value decomposition (SVD) Il27ll . As fully reconstructed jets in 
heavy ion collisions are a relatively new experimental tool, special attention is taken to discuss potential problems 
with these measurements and the measures taken to address these issues in the analysis. A measurement of the 
suppression of single inclusive muons is also presented. A template technique was used to perform a statistical 
separation of the prompt muon signal fraction from the total measured inclusive muon spectrum which includes 
experimental background from decays in flight of pions and kaons as well as muons produced due to interactions in 
the detector. 



2. Expermental setup and data sample 

The ATLAS detector consists of three main systems constructed in a barrel/end-cap geometrical confguration [28 1. 
The inner detector (ID) provides charged particle tracking over \q\ < 2.5. It is composed of silicon pixel detectors in 
the innermost layers, followed by silicon microstrip detectors and a straw-tube tracker, all immersed in a 2 T axial 
magnetic field provided by a solenoid. The ATLAS calorimeter system is composed of both electromagnetic (EM) 
and hadronic sampling calorimeters providing multiple longitudal sampling layers. These include the EM barrel and 
end-cap liquid Argon (LAr) calorimeters Qtj\ < 3.2), steel-scintillator hadronic tile calorimeter (|^| < 1.7), a LAr 
hadronic end-cap (1.5 < I77I < 3.2). The range is extended by the LAr forward calorimeters (FCal) 3.2 < \r]\ < 4.9 
which provides a combination of electromagnetic and hadronic energy measurements. Beyond the calorimeter is a 
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muon spectrometer (MS) which tracks muons in the range \rj\ < 2.7, and is composed of additional tracking chambers 
in a magnetic field produced by three air-core toroid systems. The Zero Degree Calorimeters (ZDC), positioned at 
I77I > 8.3, provide measurements of spectator neutrons and were used to define a trigger requiring a pulse height above 
a threshold set below the single neutron peak. The Minimum Bias Trigger Scintillators (MBTS) provide measurements 
of charged particles in the range 2.1 < I77I < 3.9 and consists of two pairs of counters placed at z = +3.6 m. 



The measurements presented here were performed using data from the 2010 a/^nn = 2.76 TeV Pb+Pb run cor- 
responding to a total integrated luminosity of 7pb~ l . A minimum bias sample was selected by requiring a ZDC 
coincidence trigger. In addition to the trigger requirements events were required to have a good reconstructed primary 
vertex and a measured time difference between the two sides of the MBTS of Af MBT s < 3 ns. The combination of 
trigger and offline requirements select minimum-bias hadronic Pb+Pb collisions with an efficiency of 98 + 2%. The 
resulting data sample contained approximately 50 million events after applying these requirements. The centrality 
determination was performed by measuring the minimum bias total FCal Ej, ££^ Cal , distribution and partitioning the 
range into percentile bins of the total integral. The event selection criteria and centrality definitions correspond to the 
standard ATLAS defintions for the analysis of 2010 data l29l . For both the jet and muon suppression measurements 
the 60-80% centrality bin was used as the peripheral reference in Eq. [T] 

Two separate Monte Carlo (MC) samples were used in assessing reconstruction performance and extracting cor- 
rections. The MC overlay sample was constructed by embedding PYTHIA [ 30 1 dijet events within minimum bias 
HIJING events 1311 . The separate samples of PYTHIA events were generated at different ranges of the pj in the 
hard-scattering process to insure good statistics over the full range of jet pj. These samples were used for both the jet 
and muon suppression analyses. A separate sample of minimum bias HIJING events with no embedding of PYTHIA 
jets, the HIJING only, sample was produced in an analogous fashion for the purposes of studying the absolute fake 
rate and UE fluctuations. 

3. Jet reconstruction 

The jets used in this analysis were reconstructed using the anti-£, algorithm with R = 0.2, 0.3, 0.4 and 0.5. 
The effects of the UE were removed by applying a background subtraction procedure. Calorimeter towers with 
A77 x A</> ^0.1x0.1 were constructed from calorimeter cells and then clustered using the anti-A:,. The background was 
then subtracted on a per cell basis. For the i th cell in a jet 



The average density, p, and flow modulation, v%, are determined separately for each calorimeter sampling layer j, and 
an iterative procedure is used to ensure the jets themselves do not bias the estimates of p and V2. The event plane angle 
¥2 is determined from the azimuthal Ej distribution in the FCal, which is in a separate 77 interval (3.2 < \r/\ < 4.9) 
from the jets considered here. Each cell is considered to be a four-vector using the E™. , rji and 0, with the mass taken 
to be zero. The jet's four-vector is then taken to be the four-vector sum of each of its constituent calorimeter cells. 

The energy density of the UE has contributions from the mean background as well as global correlations and local 
uncorrelated fluctuations. The subtraction procedure above removes the contribution of the mean and also the most 
significant global correlation, the elliptic flow. However, uncorrelated Ej fluctuations inside the nominal jet cone 
may significantly distort the measured jet energy on a jet-by-jet basis. This effect, combined with fluctuations in the 
energy deposition result in a jet energy resolution (JER) causing bin migration in the measured quantity. As the jet 
spectrum is steeply falling, the effect of the JER is to cause a net shift of the spectrum to higher pj. This effect must be 
corrected for by an unfolding procedure, which requires a detailed knowledge of both the detector response and the UE 
fluctuations. The jet Rqp as a function of pj for the R = 0.2 (red diamond) and R = 0.4 (black circle) jets is shown on 
the left panel of Fig.[T]both before (hollow markers) and after (solid markers) a correction for unfolding and efficiency, 
with the ratio shown below. As the UE fluctuations grow with area, the result for the smaller R = 0.2 jets changes little 
after the unfolding correction except at the lowest pj. However, the difference for the R = 0.4 jets is much greater 
and is nearly the magnitude of the measured signal, which emphasizes the requirement of having experimental control 
over the fluctuations to perform a meaningful measurement. Although the most technically challenging aspect of the 
measurement, the specifics of the unfolding procedure will not be discussed further in these Proceedings. For specific 
details of this procedure the reader is referred to Refs. [26 , [321. 
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Figure 1 : Left: measured and unfolded Sep values for the 0-10% centrality bin as a function of jet pj for S = 0.4 and R = 0.2 jets with the ratio 
of unfolded to measured Sep values for both jet radii shown below. The error bars on the points represent statistical uncertainties only. Right: 
Unfolded Sep values as a function of jet pj for R = 0.4 anti-&, jets in four bins of collision centrality. The error bars indicate statistical errors from 
the unfolding, the shaded boxes indicate unfolding regularization systematic errors that are partially correlated between points. The solid lines 
indicate systematic errors that are fully correlated between all points. The horizontal width of the systematic error band is chosen for presentation 
purposes only. Dotted lines indicate Sep = 0-5, and the dashed lines on the top panels indicate Sep = 1 1261 . 



4. Jet performance validation 

A separate study was performed to provide a quantitative comparison between the UE fluctuations in the data and 
those present in HIJING, which was used to determine the jet performance and extract the response matrices used in 
the unfolding. In this fluctuations study the HIJING only MC sample was used. The total transverse energy, "EE T , in 
rectangular "windows" of adjacent calorimeter towers was computed for all windows on the range I77I < 2.8. The study 
was performed for different window sizes, where the sizes were chosen so correspond to the nominal area of anti-£, 
jets (e.g. 3x4oS = 0.2 and 7x7oS = 0.4). The distribution of "LEj was measured in fine bins of 2£'^ Cal , and is 
shown in Fig.|2]for the 3x4, E T x4 , and 7x7, £ T 7 , window sizes in the top left and right panels respectively for both 
the data and the HIJING MC sample on the range 3.4 < ££^ Cal < 3.5 TeV, which lies in the 0-1% centrality interval. 
In both cases the mean obtained from the data has been subtracted. The results indicate that shape of the fluctuations 
spectrum produced by HIJING and that present in the data agree extremely well over several orders of magnitude. 
The tail at the high end can be attributed to jets, which were not removed from the sample. The magnitude of the UE 
fluctuations that may be present underneath a jet after the average background is subtracted is given by the width of 
this distribution. For each bin in "EE^ 1 the standard deviation of these distributions was computed for both data and 
the MC sample, with the results shown in the bottom panel of Fig. [2] This result indicates that the UE fluctuations 
present in HIJING, which are used to provide the response matrices for the unfolding procedure, are consistent with 
the data over the full centrality range. The data-MC differences obtained from the fluctuations analysis were included 
in the systematic uncertainties in the jet spectrum measurement. 

Large UE fluctuations may be falsely identified as jets, and the contribution of such "fake" jets must be understood 
and removed from the spectrum. This was achieved by requiring an additional signature of hard particle production 
in association with a reconstructed jet. For these purposes track jets, jets reconstructed from charged particles with 
pj > 4 GeV using the anti-£, algorithm with R = 0.4, and electromagnetic clusters were used. Jets were required to be 
within AR < 0.2 of either a track jet or electromagnetic cluster with pj > 7 GeV. This procedure was found to reduce 
the absolute rate of fake jets by a factor of approximately 50, while introducing only a slight additional -dependent 
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Figure 2: Top: Representative distributions of E^ x4 - <£| x4 > (left) and E\ xl - (E 1 * 1 ) (right) for data (points) and MC (filled histogram) for Pb+Pb 
collisions with 3.4 < S£^ Cal < 3.5 TeV. The vertical lines indicate E^ x4 - (Ej X4 ) = and E lxl - {E lxl ) = 0. Bottom: Standard deviations of the 
E^ xA and E lxl distributions, cr[E^ x4 ] and o-[E^ xl ], respectively, in data and HIJING MC sample as a function of S£^ Cal (26). 



inefficiency which was subsumed into the correction for the overall reconstruction inefficiency. The residual fake rate 
after rejection was estimated to be 3% of the signal rate for p T > 50 GeV. 



5. Jet suppression results 

The unfolded and efficiency-corrected jet spectra were used to construct the Rep ratios defined in Eq. [T]in bins 
of jet pi, R and centrality. The Rep for R = 0.4 jets as a function of pj is shown on the right side of Fig. [T] for 
four centrality bins. The black error bands indicate systematic uncertainties that are fully correlated among different 
Pt bins, while the shaded red boxes contain systematic uncertainties that are uncorrelated; the error bars indicate 
statistical uncertainties which develop a significant covariance through the unfolding. In the most central collisions, 
the Rep reaches a value of approximately 0.5 and shows a weak pj dependence. This pj dependence is common 
to all centrality bins, with the overall level of suppression diminishing in increasingly peripheral collisions. The R 
dependence of the suppression is shown in the left panel of Fig. [3] where the Rep ratio measured with a particular 
R value, /?* p , is divided by the corresponding value for R = 0.2 jets. As many of the uncertainties are correlated 
between the different R values, many cancel in the ratio. The statistical error bars shown account for the full statistical 
correlation between jets of different R values propagated through the unfolding procedure. The results in that figure 
indicate a significant dependence of Rep on jet radius collisions. For pi < 100 GeV the R^, p /R^p values for both 
R = 0.4 and R = 0.5 differ from one beyond the statistical and systematic uncertainties. This deviation persists for 
R = 0.5 above 100 GeV. A similar, but weaker dependence is observed in the 10-20% centrality bin, and in more 
peripheral bins, no significant radial dependence is observed. The differences between Rqp values for the different jet 
radii increase with decreasing pj, except for the lowest two pj bins. However, direct comparisons of Rqp between 
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Figure 3: Left: ratios of Rqp values between R = 0.3, 0.4 and 0.5 jets and R = 0.2 jets as a function of pj in the 0-10% centrality bin. The error 
bars show statistical uncertainties. The shaded boxes indicate partially correlated systematic errors. The lines indicate systematic errors that are 
fully correlated between different px bins | 26 1. Right: muon Sep as a function of pi in four centrality bins. The shaded boxes indicate systematic 
uncertainties that are fully correlated between pj bins. The error bars indicate combined systematic and statistical errors, which are taken to be 
uncorrelated [33). 



different jet radii at low pj should be treated with care as the same jets measured using smaller radii will tend to 
appear in lower pj bins than when measured with a larger radius. 

6. Inclusive muon analysis 

In the interval 4 < pr < 14 GeV, the inclusive production of single muon is dominated by contributions from the 
semi-leptonic decays of hadrons containing open heavy flavor. Contributions from electroweak boson decays as well 
as Drell-Yan production and quarkonia decays all contribute at less than the 1% level over this momentum range [25 1- 
Besides the heavy flavor signal, the inclusive muon yield contains backgrounds from ji and K decays inside the 
detector as well as muons produced in hadronic showers in the calorimeter. These muons have different reconstruction 
properties, which were used to construct a discriminant. The templates were constructed from composite discriminant 
distributions, dP/dC, of both signal and background sources in an MC sample and the signal fraction was extracted 
using a statistical separation. A more detailed presentation of this analysis is presented in Ref. 1551 . 

Muons are reconstructed by combining information from the ID and MS as well as a parameterized form of the 
typical muon energy loss inside the calorimeter. For background muons, which are produced within the detector, the 
momentum measurements in the ID and MS will differ by more than the typical energy loss. Therefore the fractional 
energy loss distribution, 

PlD PlD 

will differ for signal and background. Furthermore, tracks corresponding to muons produced from decays in the ID 
may exhibit kinks at the decay vertex which will exceed the typical angular deflection due to multiple scattering, </> msc . 
This can be quantified by considering the relative deflection at the k th scattering center, 



S(k) = ^= 
yn 



( k 



U=l j=k+i 



(4) 



with positions of the scattering centers coinciding with the locations of the ID hits from which the track is recon- 
structed. This was used to define a scattering significance for each track as the maximum of the individual scatterings, 
S = max [\S (k)\, k - 1,2, •■•«}. These two variables were used to construct a composite variable, 

C = E^l + rS , (5) 

PlD 
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where r is a parameter that controls the relative contribution of each variable and was fixed at r — 0.07 for optimal 
signal separation. The signal fraction, /s, was obtained by fitting the composite distribution obtained from the data, 
dP/dC, in each centrality and muon pj bin to the form 



dP dP 
dC ~* S dC 



(6) 



where df/dC| s and dP/dC| B are the composite distributions obtained from the MC for the signal and background 
contributions respectively, and were obtained by constructing template probability distribution functions using a kernel 
estimation method included in the RooFit package P4l . Template fits are shown in Fig. [4]for the 5 < pr < 6 GeV 
(top) and 10 < px < 14 GeV (bottom) in the 0-10% (left) and 60-80% (right) centrality bins. In this analysis, muons 
were required to be on the interval < 1.05. The extracted signal fraction was used to obtain the signal muon 
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Figure 4: Template fits for the composite discriminant distributions — , with C defined in Eq. 5 for muons with 5 < pj < 6 GeV (top) and 

dC □ 

10 < pi < 14 GeV (bottom) in the 0-10% (left) and 60-80% (right) centrality bins. The black points indicate the data, with the black curve 
indicating the best fit to the two-component template. The signal and background contributions are indicated by the brown filled and blue dashed 
blue curves respectively 1331 . 



spectrum from the total muon yield, yV muon , and ratios were constructed to form the Rep, 
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where s is the pj- and centrality-dependent muon reconstruction efficiency obtained from the MC sample. The 
resulting Rep values are shown in the right panel of Fig. [3] The Rep values vary weakly with p T and the points for 
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each centrality interval are consistent with a /^-independent R CP within the uncertainties on the points. The Rqp varies 
strongly with centrality, increasing from approximately 0.4 in the 0-10% centrality bin to approximately 0.85 in the 
40-60% bin. 

7. Conclusion 

Measurements of the centrality and p T dependence of the suppression of single inclusive jets have been presented 
for jets reconstructed with the anti-A;, algorithm with R = 0.2, 0.3, 0.4 and 0.5. The results show that jets are sup- 
pressed in central collisions relative to peripheral collisions by approximately a factor of two. This suppression shows 
little variation with pi over the range 38 < p? < 210 GeV for all R values. The effects of the suppression diminish in 
a gradual fashion in increasingly peripheral collisions. For jets with pj < 100 GeV, greater suppression is observed 
for jets with smaller R values at the same p T . Furthermore, the relative difference in suppression appears to be largest 
at smaller jet pj which is qualitatively consistent with expected broadening due to radiative energy loss [35]. The 
suppression of single muons from open heavy flavored hadrons was also measured and displays similar features in 
both the overall magnitude, centrality and pj dependence of the suppression. This is different than the suppression 
observed for single hadrons at comparable pj at the LHC which is rising steeply ifTOlfTTI . This feature is qualitatively 
different than the results from RHIC, where measurements of heavy quark suppression at RHIC via semi-leptonic de- 
cay channels in Au+Au collisions where comparable to the observed suppression of inclusive hadron for pj > 4 GeV 

03 nana. 
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